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ABSTRACT 

Aims. To determine whether "vortex creep" in neutron stars, the slow motion of neutron vortices with respect to pinning sites in the core or 
inner crust, is consistent with observations of long-period precession. 

Methods. Using the concept of vortex drag, I discuss the precession dynamics of a star with imperfectly-pinned (i.e., "creeping") vortices. 
Results. The precession frequency is far too high to be consistent with observations, indicating that the standard picture of the outer core 
(superfluid neutrons in co-existence with type II, superconducting protons) should be reconsidered. There is a slow precession mode, but it 
is highly over-damped and cannot complete even a single cycle. Moreover, the vortices of the inner crust must be able to move with little 
dissipation with respect to the solid. 
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There is mounting evidence that some isolated neutron stars 
undergo long-period precession (nutation). The strongest evi- 
dence is found in PSR 1828-11, which shows highly-periodic 
variations in pulse phase over a period of ~ 500 d, accompa- 
nied by coiTelated changes in beam width fetairs et al] l2000). 
PSR 1642-03 also shows periodic changes in pulse phase, 
though correlated chan ges in the beam have not been detected 
dShabanov etalJ200ll) . RX J0720.9-3125 is the first x-ray pul- 
sar to show evidence for precession; the precession period is 
~ 7 yr, with correlated changes in line depth ( Haberl et al. 
l2006l) . Low-level timing "noise", seen in all pulsars, is quasi- 
periodic in many cases, and could represent precession at 
low amplitude (for ex am ples se e, e.g. Downs & Reichlev 1983 
and Id' Alessandroet^ Il993l) . Physically-motivated models 
of precession provide good fits to the data of PSR 1828-11 
(Link &EDStein 2001; Akaunetal. 2006) and RX J0720.9- 
3125 faaberl et alJl2006l) . supporting a precession interpreta- 
tion. 

The manner in which a neutron star precesses depends on 
the dynamics of its interior, and so observations of preces- 
sion can be used to study the stellar interior The purpose of 
this Note is to discuss the role that quantized neutron vortices, 
which are expected to occupy most of the neutron star inte- 
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rior, play in the dynamics of neutron star precession. I focus 
on recent work on this question and emphasize that pinning of 
vortices to magnetic flux tubes in the core, or to nuclei in the 
inner crust, is incompatible with observations of long-period 
precession. 

The first part of this Note is aimed at constraining the vortex 
dynamics and state of the outer core. To make this argument, I 
assume a priori that there is no pinning or significant dissipa- 
tion associated with the dynamics of vortices in the inner crust. 
At the end of the Note, I argue that the assumption of weak 
dissipation in the inner crust is required by observations, and 
obtain an upper bound on how strong the dissipation can be. 

Nucleon pairing calculations predict that the outer core 
of a neutron star consists of a neutron superfluid in co- 
existence with superconduc ting protons (for a review, see, 
IPean & Hiorth-Jensenll2003|). The sup erconductor is expected 
to be in a type II state dBavm et alJ [l969a). so that the mag- 
netic field that penetrates the core is organized in flux tubes 
teavm et al.. . 1969b). long structures of microscopic cross sec- 
tion. The flux tubes probably have a very complicated arrange- 
ment that froze to the core medium when it became a supercon- 
ductor short ly after the star's formation ( Ruderman et al.,. 19981 
|jonesl2006l) . The superconducting protons, other charges {e.g., 
electrons and muons) and the crust are all coupled together 
through magnetic stresses over ti me scales of se veral seconds, 
nearly corotating as a rigid body (lEassonlll979l) . By contrast, 



2 



B. Link: Incompatibility of long-period neutron star precession with creeping neutron vortices (RN) 



the neutron superfluid rotates by establishing a nearly rectilin- 
ear array of vortex lines, whose arrangement determines the 
angular momentum of the neutron fluid. The neutron vortices 
of the outer core, w hich are themselve s magnetized through 
Fermi liquid efifects dAlpar et alJll984al) . pin to the flux tubes 
llSaulsll989l) . The origin of the pinning is that bringing a (mag- 
netized) vortex close to a flux tube raises (or lowers, depending 
on orientation) the magnetic energy by ~ 5 MeV per intersec- 
tion. For typical neutron star rotation rates and magnetic fields, 
there are ^ 10'^ vortex lines and ~ 10^' flux tubes. The vor- 
tices are thus tangled in the flux tube array, with an intersection 
spacing of ~ 10"'° cm along a vortex. This entanglement pre- 
vents the neutron fluid from corotating with the rest of the star. 

From the standpoint of precession dynamics, the star can be 
regarded as consisting of two components: 1) the core protons, 
other charges and crust, which rotate as a single "body", and, 2) 
a neutron superfluid which rotates according to the distribution 
and movement of its vortices. The angular momentum vector 
of the superfluid (most of the star) is fixed to the body to the 
extent that the vortices remain pinned against flux tubes, while, 
in a precessing body, the angular velocity vector changes its 
orientation with respect to the body's symmetry axes. As I ar- 
gue in this Note, this immobilization of vortices with respect to 
the body gives a precession period that is far to fast to reconcile 
with observations. 

Precession is a rotational mode of a rigid body in which 
the body's angular velocity is not aligned with any principal 
axis. Recall that for a rigid, biaxial body of oblateness e, the 
precession frequency m, in units of the spin frequency u, is 
vj - e. If a neutron star precesses approximately as a rigid 
body, the observed precession period of PSR 1828-11 implies 
e ^ 10"^; the other precession candidates imply similar values 
of e. 

Pinning of vortices anywhere dramatically increases the 
precession frequency. Shaham (1977) showed that if vortex 
pinning is perfectly eff'ective then precession would occur at 
very high frequency: 



m — e + 



(1) 



where Lp is the angular momentum contained in the pinned su- 
perfluid and lb is the moment of inertia of the "body" (crust 
plus charges). The pinned fluid effectively acts as a huge ad- 
ditional contribution to the star's oblateness, since Lp ^ IpU) 
and Lp/IhCij ^ Ip/h » e- The precession frequency is thus 
determined by how much of the neutron vorticity is pinned. If 
most of the neutron superfluid of the outer core is pinned to the 
flux tubes, the precession frequency will be tir ^ Ip/h - 10, 
compared to 10"^ observed for PSR 1828-11. Hence, perfect 
pinning in even a small fraction of the core gives a precession 
frequency that is far too fast to be consistent with observations. 
This result is a consequence of the fact that the vortex aiTay car- 
ries enormous vorticity; with the added constraint that there is 
vorticity fixed to the body, the star must precess at a frequency 
higher than e in order to conserve angular momentum. For pre- 
cession to occur with long period, the vortices must be able to 
closely follow the instantaneous spin vector of the rest of the 
star, and so cannot be perfectly pinned in any sizeable region 



of the star; the picture of vortex lines entangled in flux tubes 
appears to be i ncompatible with observations of long-period 
precession jLink 2003i). 

The assumption of perfect pinning is, however, an ide- 
alization. Some vortex motion with respect to pinning sites 



could occur by a process of creep through thermal activation 


(Alcar et al. 


1984b; Linketal. 1993; 


Chau & Chens 


1993) or 



quantum tunneling (ILink et al.lll993l) . lAlpail ( l2005l) has sug- 
gested that vortex creep could resolve the problem that preces- 
sion is too fast with perfectly-pinned vortices. I now show that 
even if core vortices creep with respect to the flux tubes against 
which they are pinned, the precession frequency is still too fast, 
by a factor of ~ 10^, to explain the observations. 

The problem at hand is how the creep of vortices with re- 
spect to defects affects precession. The defects I will first con- 
sider are flux tube/vortex junctions in the outer core. I will 
then consider the nuclei of the inner crust. The formulation of 
Sedrakian, Wasserman & Cordes (1999; hereafter SWC), who 
studied precession of a star with neutron vortices using a de- 
scription of vortex motion based on the concept of vortex drag, 
is particularly useful for addressing this question. They intro- 
duced a drag force per unit length of vortex of the form (SWC, 
eq. 29, with notational changes) 



fd = -'7(v,. - \b) - r]'ny X (v, - v^). 



(2) 



where v,, is the velocity of a vortex line, v;, is the velocity of the 
background body against which the vortices are dragged, 77 and 
T]' are drag coefficients, and fiy is a unit vector in the direction 
of the angular velocity of the superfluid (coincident with the 
vorticity axis and also parallel to the vortex aiTay). In the outer 
core, the background against which the vortices are dragged 
is the fence of flux tubes. The first term refers to (dissipative) 
drag anti-parallel to the vortex motion, while the second term 
describes a (possible) non-dissipative force transverse to the 
vortex motion. It is convenient to define the dimensionless drag 
coefficients 



(psK - ri'Y + if 



pMPsK - rj') 



(3) 



(4) 



where p, is the superfluid mass density and k (= h/2m„; m„ is 
the neutron mass) is the quantum of vorticity in a neutron su- 
perfluid. The drag coefficients determine how the vortex moves 
with respect to the background (SWC, eq. 33): 



V,. -\b = (J3' - l){\b - V,) + ySMv X (Vi - V,), 



(5) 



where Vj is the velocity of the neutron superfluid. 

Vortex creep, by definition, means that v,i - v;,, which 
from eqs. implies that rj/psK » 1. Hence, vortex 

creep is in the high-drag limit. To calculate the precession dy- 
namics, however, we must know the relative magnitudes of 77 
and 77'. The transverse force (second term of eq. |2} has not 
been shown to exist in any calculation of drag on vortices by 
any process in a neutron star. The existence of this term is, 
in fact, a controversial issue in laboratory superfluids. On the 
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one hand, IXhouless et alJ ( 1 19961) and IWexleil ( Il997h have ar- 
gued that this term is not present at all for an isolated vortex 
moving in an infinite, uniform superfluid (but see the com- 
ments by Hall & Hook 1998 and Sonin 1998 and the repUes 
therein). On the other hand, mutual friction experiments with 
superfluid He (for which the system is, of course, neither infi- 
nite nor uniform) do demonstrate the existence of the transverse 
force (e.g- bevan et al.Hl997i) . Whether a transverse force ex- 
ists for a vortex moving with respect to defects in a neutron 
star is unknown, so I will consider the possibility that it exists to 
maintain generality and appeal to laboratory experiments to ob- 
tain a bound. All mutual friction experiments in superfluid He 
that measure rj' find that this coefficient goes to zero for tem- 
peratures well bel ow the superfluid critical temperature {e.g., 
iBevan et"al]ll997l) . the appropriate temperature regime for a 
neutron star. As the temperature is increased, rj' jpsK usually 
increases, but never exceeds unity. Unless tj' has very different 
behavior in a neutron star, then t] » rj' in the limit of high 
drag. I henceforth assume this, but will return to the unlikely 
possibility of 77' >> 77 at the end of this Note. Eqs. Q and @ 
become 



(6) 



If there is no drag {rj - fi - fi' - Q, with 77' = 0), then v,. - \s 
and the vortices corotate with the superfluid. When the vortices 
are dragged, Vj and v,, differ Perfect pinning corresponds to 
77 —> 00, yS = 0, /?' = 1. In the high-drag limit (rj/psK » 1), 
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1 -j6' 



(7) 



SWC considered a two-component biaxial star with vortices 
moving according to the drag force of eq. and obtained the 
modes of the system. In the limit of high drag, the limit appro- 
priate to vortex creep, the modes take simple forms. There is 
a high-frequency mode of complex frequency (SWC, eq. 81, 
with ruf = -ipp) 



(8) 



The real part gives the precession frequency, and the imagi- 
nary part gives the damping rate. This mode coiTesponds to pre- 
cession at very high frequency (Ip/Ii,)(D, as found by Shaham 
(1977), with little damping (since p « 1). But there is also a 
new slow mode, of frequency (SWC, eq. 80, with tiTj = -ipd) 



m, = -e—p^ + le—p. 



(9) 



This mode can have a very low frequency compared to to. The 
mode is absent (has zero frequency) in the case of perfect pin- 
ning; it arises only if the vortices can move with respect to the 
body. Since P « 1 in the high-drag limit, as relevant to vortex 
creep, the damping rate of the mode is faster than the oscillation 
frequency by a factor of yS ' >> 1. Hence, this mode is not os- 
cillatory at all; it is highly over-damped, and never completes 
even one cycle. The problem becomes even more serious the 
more effective the pinning is (smaller /?). Hence even if the vor- 
tices move through creep, the fast precession mode of eq. (|8j 



is the only persistent and dynamically relevant one. The effects 
of t riaxial it y do n ot change the basic picture (SWC). 

(120051) concluded that the slow mode is under- 
damped based on an assumed torque between the crust and the 
superfluid of the following linear form: 



Ndi-a 



-(Wj - 6>i)/r 



(10) 



where r is a damping time, cjs is the angular velocity of the 
superfluid and cjt is the angular velocity of the body. This drag 
torque does not correctly describe the dynamics of a dragged 
vortex aiTay; the correct form, obtained by integrating r x 
from eq. Q throughout the star, takes a more complicated form 
(see SWC, eqs. 50 and 51), and is much larger than that of eq. 
J10> . Use of eq. JlOl l greatly underestimates the strength of the 
damping that arises from dissipative vortex motion, and leads 
to the incorrect conclusion that the slow precession mode is 
under-damped. 

Provided the assumption that the vortices of the crust are 
not significantly dragged remains valid, a simple and general 
conclusion follows from this analysis: long-period precession 
demands that the vortices of the outer core move with little 
drag, that is, they do not creep, but^ow essentially everywhere. 
Hence, in th e outer core, vortices and flux tubes cannot co-exist 
llLinlJl2003l) . The possibilities that might be consistent with 
long-period precession are: 1) normal neutrons, superconduct- 
ing protons, 2) superfluid neutrons, type I protons, 3) super- 
fluid neutrons, normal protons, 4) normal neutrons and protons, 
and, 5) the core magnetic field has somehow been almost com- 
pletely expelled. The existence of proton superconductivity is 
on rather firm footing, though the type of superconductivity is 
less clear. In a type I core, the magnetic flux could be orga- 
nized in slabs or other geometries of mesoscopic dimensions, 
rather than the flux tubes found in a type II superconductor. 
It might be possible for vortices to move through a type I su- 
perconductor with sufficie ntly low drag to a llow long-period, 
under-damped precession; ISedrakianI (l2005h showed that this 
scenario works for two specific geometries for the magnetic 
field. The pairing state of neutrons is far less certain. Most pair- 
ing calculations show pairing, but at least one calculation does 
not ( Schwenk & Friman 2004). Lacking a mechanism for com- 
plete expulsion of the core field, possibilities 1) and 2) seem 
the most likely within existing uncertainties regarding nucleon 
pairing. Until the issue of neutron pairing in the outer core is 
settled by first-principles calculations, astrophysical arguments 
of the type presented here are useful for constraining the vari- 
ous possibilities for the hadronic ground states. 

I now turn to the in ner cru s t, wh er e vortices could 
AlpaJ I 977t lEpstein & BavirJ 
1997t lAvQgadro et alJ l200d 



pin to lattice nuclei 

1988 ; IPizzochero et alJ , , 

IDonati & Pizzocherol |200^ . At the relatively low densi- 
ties of the inner crust, neutron pairing is well-understood and 
neutr on superfluidity there is th eoretically well-established 
(see Dean & Hiorth-Jensenll2003l for a review). Assuming, for 
example, that almost all of the inner-crust superfluid is pinned 
to nuclei and that the rest of the star rotates as a rigid body 
over time scales less than the precession period, the precession 
frequency becomes m ^ 0.01 liShaham.1977.) . again far too 
fast. There is a slow mode if the vortices creep, but as before 
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it is highly over-damped. This result does not necessarily rule 
out vortex pinning in the inner crust; it just cannot happen 
in stars that are slowly precessing. In fact, the hydrodynamic 
forces acting on pinned vortices in a star with a "wobble 
angle" 9 between the star's symmetry axis and the angular 
momentum of ~ 3°, as inferred for PSR 1828-11, would be 
sufficient to unpin the vortices in the inner crust ( Link & Cutler 
120021) . The question then becomes one of how strong these 
unpinned vortices are dragged. For the precession to be of long 
period, the vorticity axis must be able to closely follow the 
instantaneous rotation axis of the body. To do this, the vortices 
must move at a speed 

K - Vil ^RoiO^ 10"^ cms-' (PSR 1828-11). (11) 

Here R is the stellar radius and is the wobble angle. The mo- 
tion of vortices past the nuclei with which they interact excites 
waves on th e vortices (Kelvin mode s), an inherently dissipa- 
tive process. lEpstein & Baviiil ([1992") andl .Tonesl ill 99 21) studied 
this problem in the context of explaining the short spin-up time 
scale of pulsar glitches. They found Kelvin mode excitation to 
be highly dissipative for Iv,. - v/,|$ 10^ cm s"', a characteris- 
tic velocity difference that could develop as a consequence of 
vortex pinning. These treatments were done in a perturbative 
approximation, cannot be extended below Iv,. - Vi,| ~ 10^ cm 
s~\ and so cannot be applied to the precession problem; instead 
we require a description of vortex drag at much lower veloci- 
ties, in which the vortex response to motion past a nucleus is 
highly non-linear Jones ( 1998) showed that an infinitely long 
vortex moving slowly past a single nucleus tends to be tem- 
porarily trapped by the nucleus, and argued that vortices will 
pin, making it theoretically impossible (at least, within our cur- 
rent understanding of the inner crust) for a neutron star to pre- 
cess slowly Uones«2004i) . The treatment of |jonesl ( ll998l) . how- 
ever, gives an unphysical divergence in the dissip ation rate a t 
low velocities when the rate should instead vanish ( Link 2004). 
For vortex dynamics in the inner crust to allow long-period pre- 
cession, assuming no difficulties with the core of the sort de- 
scribed above, the drag coefficient must satisfy rjIp^K « hlh 
( lLink..2004: ,Sedrakiaa.2005.) ; here is the moment of inertia 
of the inner crust superfluid and h is now the moment of iner- 
tia of the crust plus all other components tightly coupled to it 
over time scales less than the precession period. If //, comprises 
most of the star, then hi Is - 100, which allows some room for 
vortex motion at moderately strong drag, rjlp^K ~ 1 - 10 for 
example. Whether or not rf can satisfy this upper bound for 

- Vil ^ 10"^ is an interesting, unsolved problem. 

Based on the above arguments, 1 conclude that in slowly 
precessing neutron stars: 

1 . vortices do not pin essentially anywhere in the star, so that, 

2. the standard picture of the neutron star outer core in which 
vortices coexist with flux tubes should be reconsidered, 
and, 

3. inner crust vortices must move with little dissipation 
(rj/psK^ 10) with respect to the solid. 

These results are robust provided that precession is real; they 
are not compromised by the numerous uncertainties regard- 
ing the composition and dynamics of the inner core. The only 



possible loophole for vortex creep to be compatible with long- 
period precession seems to be that rj' is not only finite, but for 
some reason satisfies rj' jp^K » rj/psK » 1 everywhere that 
vortices are pinned. Only under these unlikely circumstances 
can vortex creep admit a long-period, under-damped preces- 
sion mode (SWC). 

Acknowledgements. I thank A. Sedrakian and I. Wasserman for valu- 
able discussions, and the referee, P. B. Jones, for extremely helpful 
criticisms. Much of the work described here was supported by U.S. 
NSF grant AST-0406832. 

References 

Akgun, T., Link, B., & Wasserman, I. 2006, Mon. Not. Roy. Astr. Soc, 
365, 653 

Alpar, M. A. 1977, Astrophys. J., 213, 527 
Alpar, M. A. 2005, astro-ph/0505073 

Alpar, M. A., Langer, S. A., & Sauls, J. A. 1984a, Astrophys. J., 282, 
533 

Alpar, M. A., Anderson, P W., Pines, D., & Shaham, J. 1984b, 

Astrophys. J., 278, 791 
Avogadro, P., Barranco, R, Broglia, R. A., & Vigezzi, E. 2006, 

nucl-th/0602028 
Baym, G., Pethick, C, & Pines, D. 1969a, Nature, 224, 673 
Baym, G., Pethick, C, Pines, D., & Ruderman, M. 1969b, Nature, 

224, 872 

Bevan, T. D. C, Manninen, A. J., Cook, J. B., et al. 1997, J. Low 

Temp. Phys., 109, 423 
Chau, H. F. & Cheng, K. S. 1993, Phys. Rev B, 47, 2707 
D'Alessandro, P., McCulloch, P M., King, E. A., Hamilton, P A., & 

McConnell, D. 1993, Mon. Not. Roy.Astr. Soc, 261, 883 
Dean, D. J. & Hjorth- Jensen, M. 2003, Rev. Mod. Phys., 75, 607 
Donati, P & Pizzochero, P M. 2006, Phys. Lett. B, 640, 74 
Downs, G. S. & Reichley, P 1983, Astrophys. J. Suppl., 53, 169 
Easson, I. 1979, Astrophys. J., 228, 257 
Epstein, R. I. & Baym, G. 1988, Astrophys. J., 328, 680 
Epstein, R. I. & Baym, G. 1992, Astrophys. J., 387, 276 
Haberl, R, Turolla, R., Vries, C. P D., et al. 2006, A&A, 451, L17 
Hall, H. E. & Hook, J. R. 1998, Phys. Rev. Lett., 80, 4356 
Jones, P B. 1992, Mon. Not. Roy. Astr. Soc, 257, 501 
Jones, P B. 1998, Phys. Rev. Lett., 81, 4560 
Jones, P B. 2004, Phys. Rev. Lett., 92, 149001 
Jones, P B. 2006, Mon. Not. Roy. Astr. Soc, 365, 339 
Link, B. 2003, Phys. Rev. Lett., 91, 101101 
Link, B. 2004, Phys. Rev. Lett., 92, 149002 
Link, B. & Cutler, C. 2002, Mon. Not. Roy. Astr. Soc, 336, 21 1 
Link, B. & Epstein, R. I. 2001, Astrophys. J., 556, 392 
Link, B., Epstein, R. I., & Baym, G. 1993, Astrophys. J., 403, 285 
Pizzochero, P. M., Viverit, L., & Broglia, R. A. 1997, Phys. Rev Lett., 

79, 3347 

Ruderman, M., Zhu, T, & Chen, K. 1998, Astrophys. J., 492, 267 
Sauls, J. A. 1989, in Timing Neutron Stars, ed. H. Ogelman & E. P J. 

van den Heuvel (Dordrecht: Proc. NATO ASI), 457 
Schwenk, A. & Friman, B. 2004, Phys. Rev. Lett., 92, 082501 
Sedrakian, A., Wasserman, I. & Cordes, J. M. 1999, ApJ, 524, 341 

(SWC) 

Sedrakian, A. 2005, Phys. Rev. D, 71, 083003 

Shabanov, T. V., Lyne, A. G., & Urama, J. O. 2001, Astrophys. J., 552, 
321 

Shaham, J. 1977, Astrophys. J., 214, 251 
Sonin, E. B. 1998, Phys. Rev. Lett., 81, 4276 



B. Link: Incompatibility of long-period neuti-on star precession with creeping neutton vortices (RN) 



5 



Stairs, I. H., Lyne, A. G., & Shemar, S. L. 2000, Nature, 406, 484 
Thouless, D. J., Ac, P., & Niu, Q. 1996, Phys. Rev. Lett., 76, 3758 
Wexler, C. 1997, Phys. Rev. Lett., 79, 1321 



